Absence of Ferromagnetic Instability at the Metal-Insulator Transition in Si-inversion 
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We have measured the Shubnikov-de Haas oscillations in high-mobility Si MOS structures over 
a wide range of the carrier densities n > 0.77 x lO^^cm"^. This range includes the critical density 
Uc of the metal-insulator transition for two samples studied. The periodicity of oscillations clearly 
demonstrates that the electron states remain fourfold degenerate down to and at the 2D MIT. Both 
the effective spin susceptibility x* ^-^id mass m* remain finite and show no signatures of divergency 
at the critical density for both samples studied. To test possible divergency of 'X*{n) and m*(n) at 
even lower densities, we have analyzed the data on x* (f*) s.nd m* (n) in terms of a critical dependence 
X*,m* (X (n/no — Our data suggest that x* a-nd rn* may diverge at no ^ 0.5 x lO^^cm"^ 

(rs > 12), which is significantly smaller than Uc- 
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Despite intensive experimental and theoretical efforts 
(see, e.g., Ref. ||l[ for a bibhography) , the origin of the 
apparent "metal-insulator transition in two dimensions" 
(2D MIT) remains to be the subject of ongoing discus- 
sion. This phenomenon addresses a fundamental problem 
of the ground state of strongly correlated and disordered 
electron systems. In (100) Si inversion layers, the 2D 
MIT is observed at a sample-dependent critical electron 
density ric ~ 1 x 10^^ cm~^ ||]. 

One of the important unsolved problems is a possible 
magnetic instability in spin or valley systems. The 
electron-electron interactions drive a 2D system towards 
magnetic transition; numerical calculations for the crit- 
ical value of rs @] at the instability vary from 13 to 
20 H]. According to these calculations, the ferromag- 
netic transition is likely to be of first order with a com- 
plete rather than a partial (ferrimagnctic) spin polariza- 
tion. An interesting interpretation of the parallel- 
field magnetoresistance in Si inversion layers suggested a 
ferromagnetic instability at or very close to Uc [where 
ric ~ (0.8 — 0.85) X lO^'^cm"^, which corresponds to 
rs ~ 9] . The idea of magnetic instability and its possible 
link to the 2D MIT is important and requires a careful 
examination. 

In this Letter, we report on our experimental test 
of two possible scenarios for the magnetic instabil- 
ity: (i) complete spin and/or valley polarization occurs 
spontaneously at a sample-dependent critical density of 
the MIT, and (ii) x* diverges at a universal (sample- 
independent) value of n = no. To test the first sce- 
nario, we measured the Shubnikov-de Haas (SdH) oscil- 
lations and determined the degeneracy of the electron 
system across the 2D MIT in two samples (down to 
n = 0.77 X 10"cm~^). We found that the period of os- 
cillations corresponds to the double-degenerate spin and 
valley states even in the presence of external magnetic 
field B « 0.5 T. This rules out the possibility of a com- 



plete spin/valley polarization at the 2D MIT. To test the 
second scenario, we analyzed independent measurements 
of rn*in), and found that each could be de- 

scribed by the same critical dependence {n/riQ — 1)^", if 
we impose an upper limit on the density of the instabil- 
ity, no ^ 0.5 X lO^'^cm"^ (r^ > 12), and a lower limit on 
the critical index a ^ 0.6. The density no is significantly 
lower than ric for the samples studied. 

The measurements were performed on two Si-MOS 
samples: Si6-14 (peak mobihty fiP"""^ ~ 2.2m2/Vs) and 
Si5 (4.3m^/Vs), with the critical density Uc of the appar- 
ent MIT 1.0 X 10"cm-2 and 0.77 x 10"cm-^ respec- 
tively. The MOS structures are made on a (OOl)-Si wafer 
with [100] source-drain orientation; the gate oxide thick- 
ness was « 190 nm. The density of electrons was con- 
trolled by the gate voltage Vg and determined from the 
period of SdH oscillations. An in-plane field > 0.02 T 
was applied for quenching the superconductivity of the 
Al contact pads and gate electrode. Details of the exper- 
imental technique can be found in Rcfs. [pHlO|] 

Typical SdH oscillations of the resisitivity pxx = P are 
shown in Fig. 1 a as a function of B±. Due to a high elec- 
tron mobility, oscillations were detected in fields down to 
0.14 T. To examine directly the first scenario, we focus 
on the period of the SdH oscillations, the quantity which 
is not renormalized by interactions. 

For such low densities as presented in Fig. 1, the os- 
cillations p{B±) are "shaped" mostly by the spin energy 
gaps [p|,p HTl]JT^ . Figure 1 b shows that the magnitude 
of oscillations increases with an in-plane magnetic field 
_B[| . This confirms that the ratio of the Zeeman energy 
Ez = g*i-J,BB± to the cyclotron energy hujc is within the 
interval 1/2 < Ez/tiujc < 1, in good agreement with the 
measured values of x* (n) [^ , which control the calculated 
energy spectrum (the upper inset to Fig. lb). Curve 5 
in Fig. 1 a corresponds to the density n = ric for sample 
Si6-14. The latter has been determined in the insulating 
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regime, p oc exp(— A/fc^T"), by extrapolating the density 
dependence of the activation energy A(n) to zero 
(see the inset to Fig. lb). 
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FIG. 1. a) SdH oscillations for the sample Si6-14 at six 
densities near ric, T = 0.2 K. Curves 3 — 6 are terminated 
at the onset of a large insulating peak in p b) En- 

hancement of oscillations with By . The upper inset shows 
the energy spectrum for Bx = 0.5 T, n = 1.0 x lO^^cm"^, 
g*m* /2mb — 4.35, m* = 0.5me Vertical arrows depict 
spin polarization and the direction of the corresponding 
level shift with By . The lower inset illustrates determina- 
tion of the critical density for Si6-14. Densities are given 
in units of 10^^ cm"'^. 

In order to emphasize the low-field region, and to 
clearly illustrate the SdH periodicity, we present the ex- 
perimental data normalized by the amplitude of the first 
SdH harmonic Ai{B±) |l^,p. In evaluating Ai, we used 
the values of oc g*m* and m*{n) measured in 

Ref. 1^; the Dingle temperature was adjusted to match 
damping of the measured oscillations. Figure 2 shows 
oscillations of the resistivity 6p/poAi as a function of 
the Landau level filling, v = nh/eB±. The Sp{i')/ poAi 
data (dots) in Figs. 2a-c correspond to the p{B±) data 
1, 5 and 6 in Fig. 1 a. It is important to limit the field 
range B± < IT in the analysis of the SdH oscillations 
in order to diminish the magnetic-field-induced spin po- 
larization and reentrant quantum Hall effect-to-insulator 
transitions iQ. The former limitation was violated at 
u < 5: doubling of the period for v — A in Fig. 2 a illus- 
trates lifting the spin degeneracy by the perpendicular 
field B± ~ 1.3 T. The latter limitation was violated for 
1/ < 10 in Figs. 2 b-e and may account for the large oscil- 
lation amplitude. 

For all the densities studied, including n = ric (Figs. 
2 be), the low-field oscillations Sp/pQ have the period 
Ai/ = 4 that corresponds to a double-degenerate {i.e. 
unpolarized) spin and valley system. The minima of Sp 
in Figs. 2 are located at = 6, 10, 14, 18, in contrast to 



= 4, 8, 12, 16, as observed for higher densities. This is 
in agreement with earlier results [ pd|jT^ , p^ and with the 
measured x* values The sign of oscillations changes 
due to the Zeeman factor co.s{ttEz /hujc), when Ez ex- 
ceeds huJc/2 (at rs > 6.3) §,|l|l. 



< 

o 

to 



0.5 




r =7.29 










0.0 




VA 


A/V 


\ 


/ 




0.5 




g*m72m =3.37 




— n= 


1.30 



10 12 14 16 
filling factor 



18 20 22 



2 - 
■ 



< 

o 

to 



■ i',=8.32 g*m72m =4.36 — ^n=1.00 



<o 

to 




10 12 14 
filling factor 

FIG. 2. Oscillatory component of the resistivity mea- 
sured for samples Si6-14 (panels a-c) and Si5 (panels d, 
e). The data are shown as dots, the fits as lines both 
are normalized by Ai(B±). The temperature is 0.2 K for 
traces (a - c) and 0.03 K for traces d and e. The values of 
n (in units of 10^^ cm~^), r^, and (?*m* are shown in each 
panel. 

From these data, we can estimate the spontaneous spin 
polarization Pq near Uc. The total spin polarization of 
the interacting 2D electron system is: 



Pz + Po 



X*Bj_ 
2pBn 



Po- 



(1) 



Figures 2 c and 2e show that the complete spin (or val- 
ley) polarization does not occur at least down to = 6 
{Pz « 0.3), for both B,| = and B\\ = 0.3T (Fig. lb); 
thus, the spontaneous component Pq < 1 — 0.3 = 0.7. We 
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can impose even more restrictive upper limit on the spon- 
taneous spin polarization across the 2D MIT, Pq ^ 0.18, 
by noting that no nodes of beating are seen in SdH os- 
cillations in Figs. 2 c, 2 e over the interval oi v = 5 — 20. 
For example, a larger value of Pq = 0.19 would induce a 
beating node in the SdH oscillations at = 15. For even 
lower densities, the period of the p oscillations [ pl]p4| 
demonstrates that both the spin and valley states remain 
double-degenerate across the 2D MIT. These direct data 
provide very strong evidence against a complete sponta- 
neous spin/valley polarization at the sample-dependent 
critical density for n > 0.77 x lO^^cm"^. 

Below, we explore the second scenario, a divergence of 
X* and m* at a sample- independent density uq. Indeed, 
it has been found in recent measurements that both 
quantities increase with decreasing carrier concentration 
and are sample-independent (x* - to within ±2%, and m* 
- ±4%). We examine our data on x*{n) and rn*{n) [p|Jl6| 
for sample Si6-14 for possible critical density dependence 
in the form (n/ng — 1)~". For a given no, the expo- 
nents a-)^ and am have been obtained from fitting x* (") 
and m*{n). An example with uq = 0.53 x lO^^cm"^ and 
a = 0.63 is shown in Figs. 3 a, b. The exponents and 
am versus no (hash marks indicate the error bars) are 
shown in the inset in Fig. 3a. The standard deviation of 
a(no) has such a shallow minimum that the optimal no 
value could not be determined reliably; the large uncer- 
tainty was mainly caused by unknown critical range of 
densities. Similar uncertain situation was encountered in 
the critical analysis of the m* data for sample Si6-14 (see 
Figs. 3b and 3c). 

We now include into consideration the data for sample 
Si5 which has substantially lower Uc. Three additional 
m*{n) points from Si5 are shown as diamonds in Fig. 3 c. 
We find that the critical dependence cannot fit the data 
when no is taken greater than 0.65 x lO^^cm^^; this value 
sets the upper limit for possible uq. Although only the 
m*(n) data were available for Si5, x*(") can be estimated 
independently from the lineshape and phase of the SdH 
oscillations. The SdH pattern is very sensitive to the ra- 
tio Ez/hujc'. when this ratio becomes greater than 3/2 
{g*m* > 7.89) or smaller than 1/2 {g*m* < 2.63), the 
phase of oscillations changes by tt. Theoretical curves 
in Figs. 2d and 2e show that, even before the phase 
reverses for all oscillations, the highest-field oscillation 
{v = 6) splits starting from g*m* « 7. The absence of 
such behavior in the measured 5p(B±) traces enables us 
to obtain, correspondingly, the upper and lower estimates 
for g*m* at six densities in the range from n = 0.768 to 
0.884 X lO^^cm"^. Three of these estimates are shown 
by vertical bars in Fig. 3a (three others, located between 
the shown bars, are omitted for clarity). 

Small values of no ^ ric can be certainly accommo- 
dated by the critical dependence; however, we searched 
for the upper limit on no, in order to determine how close 
no could be to ric- For this reason, we have used the up- 



per limits for x* (the top of bars in Fig. 3 a) when we 
plotted the critical dependence in Fig. 3 a. The x*(n) 
data for both Si5 and Si6-14 obey a common critical 
dependence only if we choose no ^ 0.53 x lO^^cm"^. 
This choice of no also provides the lower limit for the 
critical index > 0.63. This procedure defines the 
range of densities where the critical behavior holds: e.g., 
for no = 0.53 x lO^^cm"^, this range corresponds to 
n < 1.5 X lO^^cm^^. A similar but less restrictive conclu- 
sion follows from the critical analysis for m* (Figs. 3b, c): 
no ^ 0.65 X 10^^cm~^ and > 0.4. It is important to 
note that the upper limit on no is a factor of 1.5-2 lower 
than the critical density ric for the samples studied. 
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FIG. 3. Log-log plots of (a) the spin susceptibil- 
ity X* IXb and (b,c) the mass m* /mi, vs (n/no — 1): 
dots for sample Si6-14, diamonds and bars for Si5. The 
m*(n) data are plotted for two values of no = 0.53 and 
0.65 X lO'^^ cm~^. On panels (a): the vertical bars extend 
from the upper to lower limits for x*, as discussed in the 
text. On panels (b,c): open and closed symbols depict 
the upper and lower estimate for m* , obtained from the 
T-dependence of SdH amplitude j^,^. The dashed lines 
show the critical behavior x*,"^* oc {n/no — 1)~". The 
inset shows the critical indices and am vs no. 

For a given no value, the estimated critical indices a-,(^ 
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and ttm are close to each other (see the inset to Fig. 3a). 
This might be expected: there is no experimental indi- 
cation for a critical behavior of g*(n) therefore 
X* oc g*m* and m* should exhibit the same n- depen- 
dence (a critical behavior or otherwise). 

The conclusion on the absence of the magnetic insta- 
bility at n « ric, which we have drawn from our analy- 
sis, differs from the one suggested in Refs. There 
might be several reasons for this disagreement. Firstly, 
we believe that the magnetoresistance data for in-plane 
fields, analyzed in Refs. might be indirectly related 
to the spin susceptibility of mobile electrons. Secondly, 
measurements in Refs. |^J^ were taken in strong in-plane 
fields. The strong fields drive a 2D system into the hop- 
ping regime [0; the characteristic values of go to 
zero as n approaches ric Moreover, even moderate 

fields < Ep/g* Hb induce non-linearity of magnetiza- 
tion, i.e. the x*(i?)-dependence [p|j20|. In contrast, our 
direct measurements of x* have been performed in the 
low-field linear regime. 

To summarize, we measured Shubnikov-de Haas os- 
cillations in weak perpendicular fields over a wide den- 
sity range n > 0.77 x lO^^cm"^, which includes the 
sample-specific critical densities Uc of the 2D MIT for 
two different samples. It has been found that the period 
of oscillations corresponds to the fourfold degeneracy of 
spin/valley systems on both sides of the 2D MIT. Our 
results demonstrate that the apparent 2D MIT is not 
accompanied by a spontaneous complete polarization of 
spins or valleys at zero field. Moreover, the experimen- 
tal data allow us to put an upper limit Pq < 0.18 on 
the value of a possible spontaneous polarization at the 
transition. We also explored a possibility of critical be- 
havior of the renormalized spin susceptibility and the ef- 
fective mass at a sample-independent density uq. We 
found that the divergence of both x* and m* is unlikely 
for no > 0.65 x lO^^cm"^. However, it may occur at 
lower densities (significantly less than nc): e.g. for x* ~ 
at no < 0.5 x 10"cm~2 and a > 0.6. 
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